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A  Method  for  Including  Control  Effector  Interactions  in  the  Control  Allocation 

Problem 

Michael  W.  Oppenheimer  * 

David  B.  Doman  t 

Air  Force  Research  Laboratory,  WPAFB,  OH  45433-7531 


I.  Abstract 

Much  emphasis  has  been  placed  on  over- actuated  systems  for  air  vehicles.  Over- actuating  an  air  vehicle 
provides  a  certain  amount  of  redundancy  for  the  flight  control  system,  thus  potentially  allowing  for  recovery 
from  off-nominal  conditions.  Due  to  this  redundancy,  control  allocation  algorithms  are  typically  utilized  to 
compute  a  unique  solution  to  the  over-actuated  problem.  As  the  number  of  control  effectors  placed  on  a 
vehicle  increases,  the  likelihood  of  the  occurrence  of  control  effector  interactions  increases.  For  example, 
deflection  of  an  aerodynamic  surface  that  is  upstream  of  another  aerodynamic  surface  may  cause  the  forces 
and  moments  produced  by  the  downstream  effector  to  differ  from  those  produced  when  the  upstream  ef¬ 
fector  is  not  deployed.  Another  example  can  be  found  on  launch  vehicles  that  can  use  a  combination  of 
differential  throttles  and  gimballed  nozzles  for  attitude  control.  The  effectiveness  of  gimballed  nozzles  are 
clearly  influenced  by  the  engine  thrust.  The  above  are  examples  of  the  control  effector  interaction  problem. 
In  this  work,  a  method  is  devised,  which  utilizes  linear  programming  methods  in  an  iterative  framework,  to 
take  into  account  control  effector  interactions.  While  nonlinear  programming  techniques  could  be  directly 
applied  to  such  problems,  the  lack  of  convergence  guarantees  precludes  their  use  in  flight  critical  systems. 
The  use  of  linear  programming  methods  is  appealing  because  an  optimal  solution  to  each  linear  programming 
sub-problem  can  be  determined  in  a  finite  amount  of  time. 

II.  Introduction 

Conventional  aircraft  utilize  an  elevator  for  pitch  control,  ailerons  for  roll  control,  and  a  rudder  for  yaw 
control.  As  aircraft  design  has  advanced,  more  control  effectors  (some  unconventional)  have  been  placed  on 
vehicles.  In  some  cases,  certain  control  effectors  may  be  able  to  exert  significant  influence  upon  multiple 
axes,  e.g.,  ruddervators.  When  a  system  is  equipped  with  more  effectors  than  axes  to  control,  the  system 
may  be  over- actuated  which  means  that  control  redundancy  exist.  The  allocation,  blending,  or  mixing  of 
these  control  effectors  to  achieve  some  desired  objectives  defines  the  control  allocation  problem. 

Due  to  over- actuation  and  coupling  of  control  surface  effects,  it  may  not  be  straightforward  to  determine 
how  to  translate  a  moment  or  acceleration  command  from  the  flight  control  system  into  a  set  of  commands  to 
the  control  effectors.  In  addition,  rate  and  position  limits  of  the  control  surfaces  must  be  considered  in  order 
to  achieve  a  physically  realizable  solution.  It  is  also  desirable  for  a  control  allocator  to  enable  the  aircraft  to 
recover  from  off-nominal  conditions,  such  as  a  failed  control  surface,  when  physically  possible.  A  great  deal 
of  work  in  reconfigurable  control  has  been  performed  over  the  last  decade  .1-6  Many  reconfigurable  control 
systems  make  use  of  control  allocation  algorithms  to  perform  automatic  distribution  of  the  control  power 
requests  among  a  large  number  of  control  effectors,  while  ensuring  that  the  rate  and  position  limits  of  the 
actuators  are  not  violated,  and  to  potentially  allow  recovery  from  off-nominal  conditions. 

Some  of  the  simplest  control  allocation  techniques  are  explicit  ganging,  pseudo  control,  pseudo  inverse, 
and  daisy  chaining.  Unfortunately,  each  suffers  from  difficulty  in  guaranteeing  that  rate  and  position  limits 
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will  not  be  violated  and  some  can  be  difficult  to  apply  due  to  the  need  to  derive  a  control  mixing  law  a  priori. 
Another  control  allocation  method,  called  direct  allocation,7  finds  the  control  vector  that  results  in  the  best 
approximation  of  the  command  vector  in  a  given  direction.  Unconstrained  least  squares  control  allocation 
methods,  that  account  for  rate  and  position  limits,  through  the  use  of  penalty  functions,  have  also  been 
developed.8  One  of  the  first  instances  of  linear  programming  based  control  allocators  was  from  Paradiso.9,10 
In  this  work,  Paradiso  developed  a  selection  procedure  for  determining  actuator  positions  that  was  based 
on  linear  programming  and  limited  actuator  authority.  More  recently,  the  control  allocation  paradigm  has 
been  posed  as  a  constrained  optimization  problem.11  In  this  work,  the  control  allocation  problem  was  split 
into  two  sub-problems.  The  first  was  the  error  minimization  part,  which  attempts  to  find  the  control  vector, 
such  that  the  control  effector  induced  moments  or  accelerations  match  the  desired  moments  or  accelerations. 
If  multiple  solutions  exist  to  the  error  minimization  problem,  the  second  problem  attempts  to  find  a  unique 
solution  by  driving  the  control  vector  to  a  preference  vector  and  optimizing  a  secondary  objective.  The 
linear  control  allocation  problem  has  been  extended  to  an  affine  problem12  to  account  for  nonlinearities  in 
the  moment-deflection  curves.  Quadratic  programming  has  also  been  used  in  the  past.13  An  excellent  paper 
discussing  control  allocation,  by  Bodson,14  provides  a  glimpse  into  numerous  control  allocation  techniques. 

Coupling  of  the  effects  of  control  effectors  is  a  situation  that  has  not  been  addressed  in  the  past.  In  many 
instances,  the  coupling  effects  of  aerodynamic  surfaces  are  insignificant  and  may  not  require  their  inclusion 
in  the  control  allocation  paradigm.  However,  the  coupling  effects  of  engine  thrust  and  thrust  vectoring  for 
an  ascent  vehicle  typically  cannot  be  ignored.  In  this  work,  a  method  is  developed  to  include  the  coupling 
between  two  control  effectors.  The  method  utilizes  the  simplex  algorithm  so  that  guaranteed  converge  is 
achieved.  An  example  using  fourthrust  vectored  rocket  engines  is  provided. 


III.  Linear  Control  Allocation 


The  linear  control  allocation  problem  can  be  posed  as  follows:  find  the  control  vector,  8  E  Mn,  such  that 


—  d  des 


(1) 


subject  to 


^ min  A  d  A  8max 
&  —  & max 


(2) 


where  B  E  Mmxn  is  a  control  effectiveness  matrix,  the  lower  and  upper  position  limits  are  defined  by 
8 min  £  M n  and  8 max  £  Mn,  respectively,  (5  E  Mn  are  the  control  rates,  8 max  €=  M n  are  the  maximum  control 
rates,  d des  €=  M 171  are  the  desired  moments  or  accelerations  (typically  for  inner- loop  control  laws,  d^es  £=  K3), 
n  is  the  number  of  control  effectors,  and  m  is  the  number  of  axes  to  control.  Equation  2  provides  the  position 
and  rate  limits  for  the  control  effectors.  In  a  digital  computer  implementation,  the  rate  limits  are  converted 
to  effective  position  limits.  The  combined  limits  become  the  most  restrictive  of  the  rate  or  position  limits 
and  are  specified  as 

8  <8  <8  (3) 

where 


8  = 
8  = 


min  y8maxi  8  +  A t8max 
max  (6  mim  8  At8max 


) 

) 


(4) 


Here,  8  E  Mn,  8  E  Mn,  and  8,  8  are  the  most  restrictive  upper  and  lower  control  effector  limits,  respectively. 

A  necessary  condition  for  a  system  to  be  over-actuated  is  the  number  of  columns  of  B,  n,  must  be  greater 
than  the  number  of  rows  of  B,  m.  The  true  test  of  over- actuation  is  that  the  number  of  linearly  independent 
columns  of  B  be  greater  than  the  number  of  rows  of  B.  For  inner- loop  control  laws,  the  B,  or  control 
effectiveness,  matrix  typically  becomes 


"  8L  8L  8L 
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where  L,  M,  and  TV  are  the  rolling,  pitching,  and  yawing  moments,  respectively  and  (5  =  [£i,  S2,  •  •  • ,  Sn\. 

In  some  practical  cases,  it  may  be  impossible  to  satisfy  the  equality  of  Eq.  1  subject  to  the  inequality 
constraints  of  Eq.  2.  Furthermore,  on  systems  where  control  redundancy  exists,  there  is  seldom  a  unique 
solution  to  the  problem.  In  order  to  address  these  practical  issues,  a  mixed  optimization  formulation  of  the 
control  allocation  problem,  developed  by  Bodson,14  can  be  used.  The  mixed  optimization  control  allocation 
problem  is  formulated  as  a  linear  programming  problem.  The  primary  objective  is  minimization  of  the 
difference  between  B£  and  ddes  and  the  secondary  objective  is  to  drive  the  control  effectors  to  some  preferred 
position  <Sp,  subject  to  the  constraints  of  rate  and  position  limits.  Mathematically,  the  mixed  optimization 
problem  is  posed  as: 

min  ||B<5  -  +  ||W*(<5  -  Sp)  ||i  (6) 

o 

subject  to  Eq.  3.  This  can  be  converted  to  a  standard  linear  programming  (LP)  problem: 
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where  wi,  •  •  •  ,  Wn  are  the  diagonal  elements  of  W5. 


S 
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IV.  Linear  Programming  and  Control  Effector  Interactions 

Typically,  the  effect  of  one  control  effector  on  another  is  ignored  in  the  control  allocation  problem.  In 
many  applications,  these  interaction  effects  are  small  as  compared  to  the  forces  and  moments  generated 
by  each  control  effector  acting  individually.  However,  cases  exist  where  the  interactions  should  not  be 
ignored.  In  aircraft  applications  where  the  aerodynamic  surfaces  are  used  to  modulate  forces  and  moments, 
instances  where  control  surfaces  lie  in  close  proximity  to  other  surfaces  or  cases  where  one  control  surface  lies 
downstream  of  another,  are  examples  of  situations  where  the  effectiveness  of  some  surfaces  are  influenced  by 
the  deflection  of  other  surfaces.  In  automotive  applications,  there  exists  significant  coupling  between  the  yaw 
and  axial  acceleration  effectiveness  of  the  steering  angle  and  front  wheel  braking  forces.15  Another  example 
can  be  found  on  launch  vehicles  that  can  use  a  combination  of  differential  throttles  and  gimballed  nozzles  for 
attitude  control.  The  effectiveness  of  gimballed  nozzles  are  clearly  influenced  by  the  engine  thrust.  Hence, 
a  control  allocation  method,  which  can  include  the  coupling  effects  of  multiple  control  effectors,  is  desired. 
For  an  example,  consider  the  vehicle  shown  in  Fig.  1.  This  launch  vehicle  has  control  surfaces  which  interact 
with  each  other.  In  particular,  the  flaperons  and  ruddervators.  Hence,  in  addition  to  the  conventional  forces 
and  moments  from  the  individual  deflection  of  these  surfaces,  there  are  also  forces  and  moments  from  the 
combined  deflection.  For  example,  the  pitching  moment  coefficient  is  the  sum  of  the  base  or  wing-body 
pitching  moment,  plus  the  sum  of  the  incremental  pitching  moments  produced  by  each  control  effector  taken 
one  at  a  time,  plus  the  incremental  pitching  moments  caused  by  any  interactions  between  control  surfaces: 


Cm  =  CmBase(ai  +  A  CmRF(oL,(3,M,5RF)  +  ACmiF(o;,/l,M,^f)  +  A  CmRR(a,  P,M,5Rr)+ 

ACmLR(a,P,M,  5lr)  +  ACmRF  RR(a,f3,M,  Srf,$rr)  +  ACmLF)Lfi(o;,  /?,  M,  Slf^lr) 

where  the  subscripts  RF,  LF ,  RR,  LR  imply  right  flap,  left  flap,  right  rudder,  and  left  rudder.  The  last  two 
terms  in  Eq.  9  are  the  interactions,  that  is,  a  combined  flaperon  and  rudder  deflection  cause  forces  and 
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Figure  1.  Reusable  Launch  Vehicle. 


moments.  Interactions  do  not  fit  into  a  linear  control  allocation  scheme  since  by  definition  they  are  non- 
separable  nonlinear  functions  of  two  control  deflections.  In  a  number  of  important  cases,  the  interactions 
can  be  described  by  a  bilinearity,  i.e.,  of  the  form  Cq$.  SiSj.  In  such  cases,  one  can  pose  a 

control  allocation  problem  similar  to  that  posed  in  Eq.  6  except  that  the  control  effectiveness  matrix  B  is 
replaced  by  a  control  dependent  matrix  A(S). 


min  ||A(«J)<5  -  ddes||i  +  ||W<5(<5  -  <5p)||i 

5 


(10) 


subject  to  Eq.  3,  where 
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and  QlMmMn  are  the  contributions  to  the  rolling,  pitching,  and  yawing  moments  from  the  combined 
actuation  of  2  control  surfaces.  More  specifically,  QL  can  be  expressed  as 
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The  matrices  QM  and  are  similar  to  Eq.  12.  For  the  case  where  bilinear  interaction  terms  are  the 
only  nonlinearities  of  interest,  the  main  diagonal  terms  will  be  zero.  However,  this  form  can  also  accommo¬ 
date  cases  where  forces  or  moments  are  separable  quadratic  functions  of  individual  control  deflections,  e.g., 
individual  aileron  or  flap  contributions  to  yawing  moment  at  low  angles-of- attack. 16 
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Because  the  matrix  A  itself  is  a  function  of  £,  the  control  allocation  problem  is,  strictly  speaking,  non¬ 
linear.  Rather  than  directly  applying  nonlinear  programming  techniques  to  this  problem,  we  propose  to 
iteratively  solve  a  series  of  linear  programming  (LP)  subproblems  (each  with  guaranteed  convergence  prop¬ 
erties)  in  order  to  progressively  improve  approximations  to  the  solution  of  the  original  nonlinear  programming 
problem.  Define  the  control  deflection  vector,  5  k,  as  the  most  recent  solution  to  an  LP  subproblem  that  is 
computed  as  part  of  the  iterative  procedure.  We  now  pose  the  following  LP  subproblem  whose  solution  yields 
an  updated  estimate,  5k+ 1,  of  the  control  deflection  vector  that  solves  the  original  nonlinear  programming 
problem: 

,  |X+i] 

0  0  ...  0  W\  W2  ...  wn  1  1  ...  1 


mm 

^sl)^s2  L 


Ssl 

&  s2 


(13) 


Subject  to: 


—  <^sl 

0 

—dS2 

0 

— 8k+i 

-5 

Sk+i 

< 

6 

A(^k)^k+l  —  ^s2 

d-des 

—  A(^k)^k+1  —  ^s2 

d-des 

Sk+t  —  Ssi 

6P 

&k- fvl  Ssl 

.  -dp  _ 

(14) 


The  LP  subproblems  are  solved  until  either  the  flight  control  system  requires  a  solution  at  the  end  of  a 
control  update  frame  or  until  the  following  convergence  criteria  is  satisfied: 


>k+l 


~  Sl 


<  tol 


(15) 


which  indicates  that  further  iterations  will  not  result  in  significant  improvements  in  the  solution  estimate. 

V.  Example 

In  this  section,  the  method  described  above  is  applied  to  a  hypothetical  vehicle  with  4  rocket  engines  that 
have  independent  thrust  settings  along  with  the  capability  to  gimbal  each  engine  in  the  pitch  axis.  Figure  2 
shows  the  engine  layout.  The  control  vector  is 


5  = 


3gi  d 92  ^g3  Sg4 


T\  T2  T3  T4 


(16) 


The  pitch  gimbal  angles,  denoted  by  5gi ,  are  limited  to  ±8°,  while  the  throttle  for  each  engine,  T$,  is  limited 
to  values  between  0  and  1  which  represents  a  fraction  of  the  maximum  thrust  Tmaa,  that  can  be  achieved  in 
any  engine  (for  these  initial  results,  Tmax  is  set  to  1).  It  is  assumed  that  the  distance  from  each  engine  to 
the  center  of  gravity  is  1  in  all  directions,  thus,  the  moment  arms  for  each  engine  are 


i*i  =  rlxi  +  nyj  +  rlzk  =  li  +  1  j  +  Ik 

1*2  =  r2J  +  r2y j  +  r2k  =  1  i  +  1  j  -  Ik 

1*3  =  7*3 J  +  r3y j  +  r3zk  =  1  i-  1  j  -  Ik 

r4  =  r4J  +  r4yj  +  r4zk  =  1  i-  1  j  +  lk 


(17) 


where  i,j,  k  are  unit  vectors  in  the  X,  Y,  and  Z  directions,  respectively.  The  moment  produced  by  the 
engines  is  given  by: 


4 

X 

i=  1 


M  =  ^  F  x  r  = 

i=  1 

(  TjTmax  SUK^r^)  i  ( Ti^vnax  COsJ^T*^  sin  8gi  ^ix^)  j  T  ( TjTrjHXdx  cos  5g.  Viy)  k 


(18) 
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Figure  2.  Reusable  Launch  Vehicle  Engine  Layout. 


The  components  of  Eq.  18  are  the  rolling,  pitching,  and  yawing  moments  produced  by  the  engine  throttles 
and  pitch  gimbal  angles.  Note  that  there  are  8  control  effectors  for  the  purposes  of  the  control  allocation 
problem.  The  interaction  effects  are  evident  here  as  the  products  of  throttle  and  gimbal  angle  are  included 
in  each  moment.  Entries  of  the  B  matrix  are 


dL 

ds9i 

dL 

dSjp- 


=  TjTr, 


■  COS  5 a,  Vi, 


dM  _  _  (rj- 1  rji 

d5p.  —  K-Li-Lmax 

—  —T 

—  -1  max 


—  Tm,a,x  sin  8 a  -  T i, 


dM 

dSjp. 


sill  6 q.  -\-  T^Tmax 

T, 


cos  5gir, 


cos  5girix)  ^t=TiT1 

=  T, 


T 

max-1  max 


Sin  5 a  Vi 

yx  L\ 


sin(59ir: 


dN_ 
d8 t- 


:  COS  dgji. 


(19) 

Then,  using  Eq.  5,  the  B  matrix  can  be  constructed.  For  this  example  problem,  the  form  of  the  QL  matrix 
is 


Ql  — 


where 


0 

0 

0 

0 

d2  L 

d6gidT1 

0 

0 

0 


0 

0 

0 

0 

0 

d2  L 

dd92dT2 

0 

0 


0 

0 

0 

0 

0 

0 

d2  L 

dS93dT3 

0 

d2L 


dSgzdTi 


0 

0 

0 

0 

0 

0 

0 

d2  L 

d59AdTA 

=  Tma: 


d2L 

>5gidL 

0 

0 

0 

0 

0 

0 

0 


0 

d2L 

8g2dL 

0 

0 

0 

0 

0 

0 


0 

0 

d2  L 

d893dT3 

0 

0 

0 

0 

0 


0 

0 

0 

d2  L 

d894dT4 

0 

0 

0 

0 


(20) 


(21) 


6  of  10 


American  Institute  of  Aeronautics  and  Astronautics 


In  other  words,  there  is  no  coupling  between  the  ith  gimbal  angle  and  the  jth  throttle  when  i  ^  j.  The  QM 
and  Qn  matrices  have  the  same  structure  as  QL  in  Eq.  20  with  their  entries  being 


d2M 

dSg.dTi 


-(- Trr 


d2N 

d69idTi 


sin  Sgiriz  +Trnax 

~  Tmax  SHI 


cos5girix) 


The  desired  command  vector,  ddes  is  set  to 


d des  =  ddes  +  0.05  *  rand 


(22) 


(23) 


where  rand  is  a  uniformly  distributed  random  vector  on  the  unit  interval  and  the  initial  value  for  d^es  is 


d-des  — 


.1 

.2 

.3 


(24) 


Two  problems  are  solved  in  order  to  make  a  comparison.  The  first  is  the  linear  programming  problem, 
where  the  objective  is  to  find  S  such  that  Eq.  1  holds  subject  to  Eq.  3.  In  this  case,  the  control  allocator 
does  not  have  any  information  about  the  control  effector  interactions.  If  the  interactions  are  large,  then 
the  linear  control  allocator  should  produce  a  control  vector  which  yields  large  moment  errors.  The  linear 
programming  problem  produces  Sunear  and  the  results  of  this  are  a  set  of  moments  given  by 


B  S 


linear 


(25) 


The  second  problem  solved  is  Eq.  10  subject  to  Eq.  3.  Here,  a  sequential  LP  is  solved  so  that  the  control 
effector  couplings  are  taken  into  account.  This  computation  produces  Squad  and  the  results  of  this  are  a  set 
of  moments  given  by 


1 

^quadQlj 

\ 

<  2 

SquadQM 

+  B  > 

< 

Squad^N 

(26) 


Figures  3,  4,  and  5  show  the  results  for  each  channel  of  the  problem.  For  instance,  Fig.  3  shows  the  first 
entry  in  d^es  along  with  the  the  first  entries  from  Eqs.  25  and  26.  Clearly,  the  strictly  linear  problem  is 
not  able  to  compute  the  correct  control  setting  to  achieve  the  desired  values  in  d^es.  On  the  other  hand, 
incorporating  the  interaction  terms  into  a  sequential  linear  programming  problem  yields  results  which  are 
essentially  equivalent  to  the  commands  (the  traces  from  Eq.  26  and  those  of  d^es  lie  on  top  of  each  other. 
In  an  effort  to  more  closely  show  how  the  sequential  linear  programming  problem  performs,  Fig.  6  shows 
the  results  from  all  channels,  but  does  not  include  the  linear  programming  results.  Here,  it  can  be  seen 
that  slight  differences  exist  between  the  sequential  linear  programming  results  and  the  commands,  but  the 
performance  is  clearly  superior  to  the  linear  programming  results. 


VI.  Conclusions 

The  work  presented  here  details  a  method  to  take  into  account  the  interactions  between  two  control 
effectors.  The  problem  is  solved  using  a  standard  linear  programming  solver,  however,  the  quadratic  terms 
are  included  and  the  solution  is  obtained  after  sequentially  calling  the  LP  solver.  Results  indicate  that  the 
solutions  are  superior  to  using  a  standard  linear  programming  problem. 
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Figure  4.  Linear  and  Quadratic  Control  Allocation  Results  and  ddes  ~  Channel  2. 
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Figure  6.  Quadratic  Control  Allocation  Results  and  ddes  ~  All  Channels. 
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